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Abstract 



We introduce a new set of generalized Fokker-Planck equations that con- 
serve energy and mass and increase a generalized entropy functional until a 
maximum entropy state is reached. Nonlinear Fokker-Planck equations asso- 
ciated with Tsallis entropies are a special case of these equations. Application 
of these results to stellar dynamics and vortex dynamics are proposed. Our 
prime result is a novel relaxation equation that should offer an easily imple- 
mentable parametrization of 2D turbulence. Usual parametrizations (includ- 
ing a single turbulent viscosity) correspond to the infinite temperature limit 
of our model. They forget a fundamental systematic drift that acts against 
diffusion as in Brownian theory. Our generalized Fokker-Planck equations can 
have applications in other fields of physics such as chemotaxis for bacterial 
populations. We propose the idea of a classification of generalized entropies 
in "classes of equivalence" and provide an aesthetic connexion between topics 
(vortices, stars, bacteries,...) which were previously disconnected. 



PACS numbers: 05.90.+m, 05.70.-a, 47.10. +g, 47.32.-y 
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I. INTRODUCTION 



The statistical mechanics of systems with long-range interactions is currently a topic 
of active research [1]. Systems with long-range interactions are numerous in nature: 
self-gravitating systems, two-dimensional vortices, non-neutral plasmas, metallic clusters, 
dipoles, fracture etc... These systems exhibit similar features such as negative specific heats, 
inequivalence of statistical ensembles, phase transitions and self-organization. Since they 
are non-extensive and non-additive, the construction of an appropriate thermodynamics is 
a challenging problem. Among all the previous examples, self-gravitating systems and 2D 
vortices play a special role because they both interact via an unshielded Newtonian potential 
(in dimensions D = 3 or D = 2) and possess a rather similar mathematical structure [2] . 

It has been recently argued that the classical Boltzmann entropy may not be correct 
for systems with long-range interactions and that Tsallis entropies, also called g-entropies, 
should be used instead [3]. In the context of 2D turbulence, Boghosian [4] has interpreted 
a result of plasma physics [5] in terms of Tsallis generalized thermodynamics. In the as- 
trophysical context, Plastino & Plastino [6] have noted that the maximization of Tsallis 
entropies leads to stellar polytropes, thereby avoiding the infinite mass problem associated 
with isothermal systems obtained by maximizing the Boltzmann entropy. However, the 
arguments advocated to justify Tsallis entropies in the context of 2D turbulence and stel- 
lar dynamics are usually unclear and misleading and were criticized in our previous papers 
[7-9]. In particular, we have argued that Tsallis entropies are particular if-functions (not 
true entropies) [10] whose maximization at fixed mass and energy determines (nonlinearly) 
dynamically stable stationary solutions of the 2D Euler or Vlasov-Poisson systems [11,12]. 
The if-functions can be useful in describing the metaequilibrium states resulting from incom- 
plete violent relaxation [13]. In this context, the aforementioned maximization problem is 
a condition of dynamical stability, not a condition of thermodynamical stability. Therefore, 
Tsallis entropies have not a fundamental justification for stellar systems and 2D turbulence. 
They form just a one-parameter family of if-functions that leads to simple models (stellar 
polytropes and polytropic vortices). Tsallis distributions can provide, however, a conve- 
nient fit of the metaequilibrium state in case of incomplete relaxation. In that context, the 
parameter q, which can vary in space, measures the importance of mixing [9]. 

On the other hand, it has been shown that Tsallis generalized thermodynamics could 
be useful to interpret anomalous diffusion in complex systems and that the g-entropies are 
connected to nonlinear Fokker-Planck equations [14]. In fact, Tsallis entropies are just a 
particular class of a much larger class of functionals that we shall call generalized entropies. 
These functionals are defined as S = — J C(f)d 3 rd 3 v where C(f) is a convex function 
of the distribution function. Many important properties obtained with Tsallis g-entropies 
(nonlinear Fokker-Planck equations, generalized if -theorem, Legendre transforms,...) re- 
main valid for these more general functionals. Tsallis entropies give a special importance 
to power-laws. Power-laws are indeed important in physics (in relation, among other, with 
multifractality [3]) but they are not the most general distributions. Generalized entropies 
arise naturally when the diffusion coefficient is an arbitrary function D(f) of the distribution 
function. When the diffusion is counter-balanced by a friction or a drift, they play the role 
of Lyapunov functionals and satisfy a iJ-theorem S > 0. In this context, Tsallis entropies 
correspond to a power-law dependance of the diffusion coefficient D(f) ~ f 9 ' 1 and the 
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g-parameter in Tsallis formalism is related to the exponent of anomalous diffusion. 

In the first part of the paper, we develop a generalized thermodynamical formalism for 
a large class of entropy functionals encompassing Boltzmann, Fermi and Tsallis entropies. 
This formalism can have applications in different domains of physics (or biology, economy, 
mathematics,...) with different interpretations. In Sec. II A, we introduce generalized 
entropies that extend those introduced by Tsallis and co-workers. In Sec. II B, we establish 
the conditions of generalized thermodynamical stablity in microcanonical and canonical 
ensembles and discuss the possible inequivalence of statistical ensembles when the caloric 
curve presents turning points or bifurcations. In Sees. II C and II D, we introduce generalized 
kinetic equations (Fokker-Planck, Landau) that conserve mass and energy and increase a 
generalized entropy functional instead of the Boltzmann entropy. In Sec. Ill, we study the 
generalized Smoluchowski-Poisson system and mention possible applications to bacterial 
populations (chemotaxis). 

In the second part of the paper (Sec. IV), we discuss the statistical mechanics of Hamil- 
tonian systems with long-range interactions focusing on stellar systems and 2D point vor- 
tices (or inviscid continuous vorticity fields). We consider the collisionless regime when the 
N — > +00 limit is taken before the t — > +00 limit (Vlasov limit). Due to mean- field effects 
and long-range interactions, the system undergoes a violent relaxation (Sec. IV A). The 
resulting metaequilibrium state is a stationary solution of the Vlasov (or Euler) equation on 
the coarse-grained scale. Its nonlinear dynamical stability can be settled via a thermody- 
namical analogy (Sec. IV B). Generalized entropies (also called H-functions) arise due to the 
existence of fine-grained constraints (Casimirs), non-ergodicity and non-ideal effects (forcing, 
dissipation,...). We propose a novel relaxation equation (Sec. IV C) that can serve either 
as a small-scale parametrization of 2D turbulence (Sec. IV E) or as a powerful numerical 
algorithm to compute arbitrary nonlinearly dynamically stable stationary solutions of the 
2D Euler- Poisson system (Sec. IV D). We propose the idea of a classification of generalized 
entropies in "classes of equivalence" with the heuristic argument that entropies of the same 
class should lead to similar results (Sec. IV F). 



Let us consider a system of N particles in interaction and denote by /(r, v, i) their 
distribution function defined such that fd 3 rd 3 v gives the total mass of particles with position 
r and velocity v at time t. Let F(r,t) = — V$ be the force (by unit of mass) experienced 
by a particle. We assume that the potential $(r, t) is related to the density p(r,t) = J fd 3 v 
by a relation of the form $(r) = / p(r')u(r — r')d 3 r' where u(r — r') is an arbitrary binary 
potential. For example, if u(r — r') = —G/\r — r'|, $ is solution of the Poisson equation 



II. GENERALIZED THERMODYNAMICS AND FOKKER-PLANCK 

EQUATIONS 



A. Generalized entropies 



A$ = 4nGp. 



(1) 



We assume that the system is isolated so that it conserves mass 




(2) 
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and energy 



E = J ^fv 2 d 3 rd 3 v + \jf& d 3 r = K + W, (3) 

where K is the kinetic energy and W the potential energy. The conservation of angular 
momentum L = / /(r x v)d 3 rd 3 v and linear impulse P = / fvd 3 rd 3 v can be easily incorpo- 
rated in the formalism. The following results remain valid if $ = $ ext (r) is a fixed external 
potential, in which case W = J p& ext d 3 r. 

We introduce a generalized entropy of the form 

S = ~J C(f)d 3 rd 3 v, (4) 

where C(f) is a convex function, i.e. C"(f) > 0. We are interested by the distribution func- 
tion / which maximizes the generalized entropy (4) at fixed mass and energy. Introducing 
appropriate Lagrange multipliers and writing the variational principle in the form 

5S - (35E - a5M = 0, (5) 

we find that the critical points of entropy at fixed mass and energy are given by 

C(f) = -j3e - a, (6) 

2 

where e = \ + $ is the energy of a particle by unit of mass. The Lagrange multipliers f3 and 
a are the generalized inverse temperature and the generalized chemical potential. Equation 
(6) can be written equivalently as 

f = F(Pe + a), (7) 

where F(x) = (C") _1 (— x). From the identity 

/'(e) = -P/C"(f), (8) 

resulting from Eq. (6), /(e) is a monotonically decreasing function of energy if f3 > 0. The 
conservation of angular momentum can be easily included in the variational principle (5) 
by introducing an appropriate Lagrange multiplier fl. Equation (6) remains valid provided 
that e is replaced by the Jacobi energy e,j = e — fi • (r x v) = |(v — ft x r) 2 + $ e jj where 
& e ff = $ — |(fi x r) 2 is the effective potential accounting for inertial forces. 

Among all functionals of the form (4), some have been discussed in detail in the literature. 
The most famous functional is the Boltzmann entropy 

S B [f] = -j f\nfd 3 rd 3 v. (9) 

It leads to the isothermal (or Boltzmann) distribution 

/ = Ae-* (10) 
Closely related to the Boltzmann entropy is the Fermi-Dirac entropy 
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Sf.d.U] = - /(-In - + fl - -) Infl - ^))d 3 rd 3 v, (11) 

which leads to the Fermi-Dirac distribution function 

f = ^ (12) 

The Fermi-Dirac distribution function (12) satisfies the constraint / < r] which is related to 
Pauli's exclusion principle in quantum mechanics. The isothermal distribution function (10) 
is recovered in the non-degenerate limit / <C i]q. We can also consider the case of bosons 
with the sign + in Eq. (11). Recently, there was a considerable interest for functionals of 
the form 

S q [f] = -^- [ J(f q -f)dVv, (13) 

where q is a real number. Such entropies introduced by Tsallis [3] are called g-entropies. 
They lead to "polytropic" distributions of the form 

f = A(X-e) n -l, (14) 
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with A = [(q — l)(3/q]~ and A = [1 — (q — l)a]/(q — l)(3. The index n of the polytrope 
is related to the parameter q by the relation n = 3/2 + l/(q — 1). Isothermal distribution 
functions are recovered in the limit q — > 1 (i.e. n — > +00). 

In any system with / = /(e), one may define a local energy dependant excitation tem- 
perature by the relation 

1 d,nf (15) 



T(e) de 

For the isothermal distribution (10), T(e) coincides with the thermodynamic temperature 
T = 1/(3. For the polytropic distribution (14), T(e) = (q — 1)(A — e). This excitation 
temperature has a constant gradient dT/de = q — 1 related to Tsallis g-parameter (or 
equivalently to the index n of the polytrope). The other parameter A is related to the value 
of energy where the temperature reaches zero. 



B. Generalized thermodynamical stability 

In the preceding section, we have just determined critical points of the generalized entropy 
(4) by cancelling its first order variations with appropriate constraints. We now turn to the 
thermodynamical stability of the solutions (in the generalized sense). We must select maxima 
of S[f] at fixed mass and energy. The condition that / is a maximum of S at fixed mass and 
energy is equivalent to the condition that 5 2 J = 5 2 S — (35 2 E is negative for all perturbations 
that conserve mass and energy to first order. This condition can be written 

5 2 J=- J C"(f)^^d 3 rd 3 v ~\(3 j 5p5^d 3 r < 0, 

V 8f I 5E = 5M = 0. (16) 
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So far, we have implicitly worked in the microcanonical ensemble in which the energy is 
fixed. However, it may be of interest to study in parallel the canonical ensemble in which 
the temperature T = 1/(5 is fixed instead of the energy. In that case, the appropriate 
thermodynamical potential is the free energy F = E — TS that we write for convenience in 
the form of Massieu function 

J = S-(3E. (17) 

According to Eq. (5), we have 

5 J = -E5{3 + a5M. (18) 

Therefore, the equilibrium state in the canonical ensemble is a maximum of J at fixed mass 
and temperature. If we just cancel the first order variations of J, this again yields the 
relation (6). The condition of thermodynamical stability in the canonical ensemble requires 
that / is a maximum of J at fixed mass and temperature. This is equivalent to the condition 
that 5 2 J is negative for all perturbations that conserve mass. This can be written 

C "(f)^£fd 3 rd 3 v - \(3 [ 5p5$d 3 r < 0, 

V5f\ 5M = 0. (19) 

We note that canonical stability implies microcanonical stability but the converse is wrong 
in general. Indeed, if inequality (19) is satisfied for all perturbations that conserve mass, it 
is a fortiori satisfied for perturbations that conserve mass and energy. Since the converse is 
wrong, this implies that we can "miss" some relevant solutions by working in the canonical 
ensemble instead of the microcanonical one. 

For self-gravitating systems described by the Boltzmann entropy (9), it is well-known 
that the statistical ensembles are non-equivalent [15,2,9]. Indeed, an isothermal distribution 
(10) can be stable in the microcanonical ensemble (maximum of Sb at fixed M and E) but 
unstable in the canonical ensemble (minimum or saddle point of Jb at fixed M and T). 
In fact, the inequivalence of statistical ensembles for systems with long-range interactions 
is not limited to self-gravitating systems nor to the Boltzmann entropy (9). It occurs for 
many other physical systems and for various functionals of the form (4). There will be 
inequivalence of statistical ensembles when the caloric curve (3(E) presents turning points 
leading to regions of negative specific heats, or said differently, when the entropy S(E) has a 
convex dip [16,17]. The stability of the solutions can be decided by using the turning point 
criterion of Katz [18] which extends the theory of Poincare on linear series of equilibria. It 
is found that a change of stability in the series of equilibria occurs in the microcanonical 
ensemble when the energy is extremum and in the canonical ensemble when the temperature 
is extremum. Stability is lost or gained depending on whether the series of equilibria turns 
clockwise or anti-clockwise at that critical point. A change of stability along a series of 
equilibria can also occur at a branching point [18,19], where the solutions bifurcate. A 
general classification of phase transitions for systems with long-range interactions has been 
proposed recently by Bouchet & Barre [20] . 
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C. Generalized Kramers equation 



We shall now introduce a generalized Fokker-Planck equation, consistent with the ther- 
modynamical framework developed previously, by using a Maximum Entropy Production 
Principle [21]. To apply the MEPP, we first write the relaxation equation for the distribu- 
tion function in the form 

| + U 6 .V 6 / = -^, (20) 

where U6 = (v, F) is a generalized velocity field in the six-dimensional phase space {r, v}, 
V6 = (d/dr, d/dv) is a generalized gradient and J/ is the diffusion current to be determined. 
The form of Eq. (20) ensures the conservation of mass provided that J / decreases sufficiently 
rapidly for large |v|. From Eqs. (3), (4) and (20), it is easy to put the time variations of 
energy and entropy in the form 

E = J J r vd 3 rd 3 v, (21) 



S = -jc"(f)J r ^-dhd 3 v, (22) 

where we have used straightforward integrations by parts. Following the MEPP, we shall 
now determine the optimal current J/ which maximizes the rate of entropy production (22) 
while satisfying the conservation of energy E = 0. For this problem to have a solution, we 
shall also impose a limitation on the current | J/|, characterized by a bound C(r, v, t) which 
exists but is not known, so that 

0<C(r,v,f). (23) 

It can be shown by a convexity argument that reaching the bound (23) is always favorable 
for increasing S, so that this constraint can be replaced by an equality. The variational 
problem can then be solved by introducing at each time t Lagrange multipliers (3 and 1/D 
for the two constraints. The condition 

SS - p(t)6E - J ^(^j)d 3 rd\ = 0, (24) 



yields an optimal current of the form 



3, = -D 



fC"(f)^ + (3(t)fv 



(25) 



The time evolution of the Lagrange multiplier (3(t) is determined by the conservation of 
energy E = 0, introducing Eq. (25) in the constraint (21). This yields 

_ S DfC»{f)%.v<PT<Pv 
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Note that the optimal current (25) can be written J/ = Dfda/dv where 



a(r,v,t) = -C'(f)-(3e, 



(27) 



is a generalized potential which is uniform at equilibrium according to Eq. (6). Therefore, 
the MEPP is just a variational formulation of the linear thermodynamics of Onsager. 

Introducing the optimal current (25) in Eq. (20), we obtain the generalized Fokker- 
Planck equation 



at 



D 



fC"(f)jt+l3(t)fv 



(28) 



Morphologically, Eq. (28) extends the usual Kramers equation introduced in the context of 
colloidal suspensions [22] and collisional stellar dynamics [23]. The first term is a generalized 
diffusion (depending on the distribution function) and the second term is a friction. The 
function j3(t) can be considered as a time dependant inverse temperature evolving with time 
so as to conserve energy (microcanonical formulation). The friction coefficient £ = Df3 
satisfies a generalized Einstein relation. Note that D is not determined by the MEPP since 
it is related to the unknown bound C(r, v, t) in Eq. (23). We can use this indetermination 
to write Eq. (28) in the alternative form 



D' 



df /3(t) 
L^v C"(f) 



(29) 



which will have the same general properties as Eq. (28). This equation involves an ordinary 
diffusion and a nonlinear friction. Equation (29) can be deduced from Eq. (28) by the 
substitution D' = DfC"(f). One of these two forms will be prefered depending on the 
situation contemplated. Note that D or D' can depend on r, v, t without altering the general 
properties of the equations. 

It is straightforward to check that Eq. (28) with the constraint (26) satisfies a if-theorem 
for the generalized entropy (4). From Eqs. (22) and (25), we can write 



S = - 



f 



d 3 rd 3 v + (3{t) J J f ■ vd 3 rd 3 v. 



(30) 



The last quantity vanishes due to the conservation of energy (21). Therefore, 



which is positive provided that D > 0. Now, at equilibrium S 
according to Eq. (25), 



0, hence J/ 



dC'(f) 



<9v 



+ /3v = 0. 



(31) 
0, so that 

(32) 



Integrating with respect to v, we get 



C(f) = -P- + A(r). 



(33) 
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The cancellation of the advective term U6 ■ V6 in Eq. (28) combined with Eq. (33) implies 
that / = /(e) and VA = — /3V$. Therefore, A(r) = — /3<&(r) — a and we recover Eq. (6) 
with (5 = lim t ^ +00 f3(t). Therefore, a stationary solution of Eq. (28) extremizes the entropy 
at fixed energy and mass. In addition, only maxima of S at fixed M and E are linearly 
stable with respect to the generalized Fokker-Planck equation (28). Indeed, considering the 
linear stability of a stationary solution of Eqs. (28) and (26), we can derive the general 
relation (see Appendix A) 

2X5 2 J = 5 2 S > 0, (34) 

connecting the growth rate A of the perturbation 5f ~ e xt to the second order variations 
of the free energy J = S — f3E and the second order variations of the rate of entropy 
production 8 2 S > 0. Since the product \5 2 J is positive, we conclude that a stationary 
solution of the generalized Fokker-Planck equation (28) is linearly stable (A < 0) if and only 
if it is an entropy maximum at fixed mass and energy (see Sec. II B). This aesthetic formula 
shows the equivalence between dynamical and thermodynamical stability for our generalized 
Fokker-Planck equations. Therefore, they only select maxima of S, not minima or saddle 
points. 

A relaxation equation appropriate to the canonical situation can be obtained by maxi- 
mizing J = S — f3E with the constraint (23). The variational principle 

SJ-J ^(0)rf 3 rrf 3 v = O, (35) 

again yields an optimal current of the form (25) but with constant (5. Since 

j=-Jj r C"{f)^j- + /3v d 3 rd 3 v = J irf 3 rrf 3 v > 0, (36) 

according to Eqs. (21), (22) and (25), we find that the free energy J increases monotonically 
until an equilibrium state of the form (6) is reached. In the canonical ensemble, we can show 
that 2\5 2 J = 5 2 J > and conclude that a stationary solution of the generalized Fokker- 
Planck equation (28) with constant j3 is linearly stable if and only if it is a maximum of free 
energy at fixed mass and temperature (see Sec. II B). 

We can also use the MEPP to construct a more general relaxation equation. Assuming 
that the diffusion current in Eq. (20) depends on v and r and repeating the same steps as 
before, we get 

df 



-1 + \J 6 -V 6 f = V 6 { j d[/C"(/)V 6 / + /3(t)/U 6± ] }, (37) 



with 



m = - ■ (38) 

where U 6 ^ = (— F,v). The generalized velocity U 6 in phase space is very similar to the 
velocity field of a two-dimensional incompressible fluid (see Sec. IV A). We note in particular 
that U 6 ^ = Vq€, where e = \ + $ plays the role of a generalized streamfunction. 
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To conclude this section, it can be of interest to discuss some special cases explicitly. For 
the Boltzmann entropy (9), C"(f ) = 1/f and Eq. (28) has the form of an ordinary Kramers 
equation 



df 



(39) 



For the Fermi-Dirac entropy (11), C"(f) = l/f(i] — /). In order to avoid the divergence of 
the term fC"(f) as / — » r]o, it is appropriate to consider the alternative form (29) of the 
generalized Kramers equation. This yields 



df 
dt 



+ U 6 • V 6 / 



dv 



+ Mvo-f)v 



(40) 



which has been initially proposed in [21]. Finally, for the Tsallis entropy (13), C"(f) = qf q 2 
and Eq. (28) has the form of a nonlinear Fokker-Planck equation 



dt 



+ U 6 • V 6 / 



d_ 
dv 



D 



dp_ 
dv 



+ /3/v 



(41) 



This equation has been studied in detail recently in relation with Tsallis entropy and anoma- 
lous diffusion [14]. The underlying mechanism giving rise to anomalous diffusion may differ 
depending on the physical system: Levy walkers, porous media, vortex dipoles in 2D turbu- 
lence, etc... In such systems, the diffusion coefficient ~ ] q ~ x is a power-law of the distribution 
function and the phase space has a fractal or multifractal structure (the exponent q is re- 
lated to the fractal dimension). In fact, the nice properties of Eq. (41), in particular the 
H-theorem, go beyond the form of entropy considered by Tsallis and remain valid for all 
convex function C(f) even if the results are not always explicit. In the context of anomalous 
diffusion, Eq. (28) can be obtained from a Langevin equation of the form 




R(t), 



(42) 



where R(t) is a white noise. Since the function in front of R(t) depends on (r, v), the last 
term in Eq. (42) is a multiplicative noise. When C(f) is a power law, Eq. (42) reduces to 
the stochastic equations studied by Borland [24]. 



D. Generalized Landau equation 



We shall now introduce another relaxation equation tending to a state of maximum 
generalized entropy at fixed mass and energy. This is the generalized Landau equation 



Sv'K^ff 



c"(f)^--c"(f) df 



dv v 



dv' v 



(43) 



u \ 



u^u v 



(44) 
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where / = /(r, v, i), f = f(r,v',t), u = v' — v, A is a constant and C(f) in any convex 
function. We shall also consider the alternative form 



| + u 6 .v 6 / = A 



/ 



1 df 1 df 



C"(f) dv» C"(f) dv' v 



(45) 



Morphologically, these equations extend the usual Landau equation introduced in plasma 
physics and collisional stellar dynamics [25,26]. The generalized Landau equation satisfies 
the conservation of mass, energy, angular momentum and linear impulse and increases a 
generalized entropy (H-theorem). In addition, formula (34) remains valid so that a stationary 
solution of the generalized Landau equation is linearly stable if and only if it is a maximum 
of the generalized entropy (4). 

Let us consider special cases explicitly. For the Boltzmann entropy, Eq. (43) reduces to 
the usual Landau equation 



dj_ 

dt 



For the Fermi-Dirac entropy, Eq. (45) takes the form 



or 



dv' v 



Finally, for the Tsallis entropy, we get 



df +U 



v 6 / = 



d 



,df x .df'« 



dv h 



f 



dv' 1 



(46) 



(47) 



(48) 



This could be called the g-Landau equation. Contrary to the nonlinear Kramers equation 
(41), it seems that Eq. (48) has never been introduced previously. We shall study its 
properties more specifically in a future work [27]. A connexion between the (generalized) 
Landau equation and the (generalized) Kramers equation can be found in a thermal bath 
approximation. In this context, / = f(v,t) describes a test particle and /' = f(y',t) 
describes the field particles. If we replace /' in Eq. (43) by its equilibrium value (6), then 
Eq. (43) becomes equivalent to the generalized Kramers equation (28) and the diffusion 
coefficient can be calculated [27]. 



III. THE GENERALIZED SMOLUCHOWSKI-POISSON SYSTEM 

A. The high friction limit 

The Kramers-Poisson system (39)-(l) is relatively complicated because it has to be solved 
in a six- dimensional phase space. However, it is well-known in Brownian theory [22] that, 
in the high friction limit £ — > +oo (or equivalently for large times t ^> £ _1 ), the velocity 
distribution function becomes close to the Maxwellian distribution and the evolution of the 
spatial density p(r, t) is governed by the Smoluchowski equation 

(49) 



dt 



-(TVp + pV$) 
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The Smoluchowski-Poisson system has been studied in [28,29]. It models the dynamics of 
self-gravitating Brownian particles and the chemotactic aggregation of bacterial populations. 

We now proceed in deriving a generalized Smoluchowski equation by taking the high 
friction limit of the generalized Kramers equation. We shall assume that (3 is constant 
(canonical situation). Integrating Eq. (28) over velocity, we get the continuity equation 

f + V( P u) = 0, (50) 

where u = (1/p) / fvd 3 v is the local velocity. Multiplying Eq. (28) by v and integrating 
over velocity, we get the momentum equation 



(p Ul ) + — ( P u lUj ) + — Pij + p — = -jD- " - ' 



dt dxj 



fC"(f)^ + pfv x 

OVi 



d 3 v, (51) 



where Pij = J fwiWjd 3 v is the stress tensor and w = v — u the relative velocity. Introducing 
the notation <f)(f) = J xC"(x)dx, the first term in the collision term can be rewritten 
d(f)(f)/dv and, since it is a gradient of a function, it vanishes by integration. We are left 
therefore with 

d d <9 5$ 

dt {fMi) + dx- {pU ^ + Wf* + P dx~ = ~ DPpUi - (52) 

We shall close this hierarchy of equations with the Local Thermodynamical Equilibrium 
(LTE) condition 

C'(f) = -p(^- + \(r,t)), (53) 

where A(r, t) depends only on position and time. This distribution function maximizes the 
local density of free energy at fixed density p and velocity u. The Lagrange multiplier A is 
related to the spatial density p(r, t) through the relation 



= / fd 3 v. (54) 



Furthermore, since the velocity distribution in Eq. (53) is isotropic, we have = pSij 
where p(r, t) is the local pressure 

p = l -J fw 2 d\, (55) 

determined by Eq. (53). From these two relations, we find that the fluid is barotropic in 
the sense that p(r, t) = p[p(r, t)] where the function p[ ] is completely specified by C(f). We 
have thus obtained what might be called the damped Euler- Jeans equations 

f + V( P u) = 0, (56) 
P< ^ = ~ VP ~ PV$ " ^° U ' (57) 
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i(Vp + pV$) 



(59) 



where d/dt — d/dt + u • V is the material derivative. These equations were first proposed in 
[21]. In the high friction limit, these equations can be simplified further since, to first order 
in we have 

pu = _I(V P + pV<£), (58) 

which is obtained from Eq. (57) by neglecting the advective term. Note that the high friction 
limit is consistent with the Local Thermodynamical Equilibrium (53). Indeed, in the high 
friction limit £ = Df3 — > +oo, the term in bracket in Eq. (28) must vanish so that the 
distribution function satisfies (53) in good approximation with u = 0(£ _1 ). Inserting the 
relation (58) in the continuity equation (56), we get a generalized form of the Smoluchowski 
equation 

dp 
dt 

initially proposed in [21]. The aesthetic form of this equation, in which the pressure p replaces 
naturally the usual term pT in the familiar Smoluchowski equation, suggests by itself the 
consistency of the generalized thermodynamical formalism. The condition of stationarity in 
Eq. (59) corresponds to an equilibrium between the pressure force — Vp and the mean-field 
force — pV^. This is equivalent to Eq. (6). Indeed, using Eqs. (54), (55) and (6), one has 

so that we obtain the condition of hydrostatic equilibrium 

Vp = -pV$. (62) 

Finally, we can show that the generalized Smoluchowski-Poisson system (59) satisfies a form 
of Virial theorem 

±Z^L = 2K + W-3p b V, (63) 

where I = J pr 2 d 3 r is the moment of inertia and pb the pressure on the box (assumed 
uniform). The proof is the same as that given in Appendix D of Ref. [29]. 

To conclude this section, we can consider particular forms of the generalized Smolu- 
chowski equation. For the Boltzmann entropy (9), Eqs. (53), (54) and (55) lead to the 
isothermal equation of state p = pT and to the usual form (49) of the Smoluchowski equa- 
tion. The case of the Fermi-Dirac entropy (11) has been treated in [21]. For the Tsallis 
entropy (13), Eqs. (53), (54) and (55) lead to the polytropic equation of state p = Kp 1 with 
7 = 1 + 1/n and Eq. (59) becomes the nonlinear Smoluchowski equation 



dt 



^(KVp 7 + pV$) 



(64) 



The nonlinear Smoluchowski-Poisson system has been studied in detail in Ref. [30]. It 
describes self-gravitating Langevin particles experiencing anomalous diffusion. It is likely 
that anomalous diffusion occurs also in biological systems so that nonlinear Smoluchowski 
equations can find applications in the context of chemotaxis. 
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B. The Lyapunov functional 

In Sec. II C, it was indicated that the Lyapunov functional associated with the general- 
ized Kramers equation (28) with a fixed inverse temperature (3 is the free energy J = S — (3E 
with the generalized entropy (4). Therefore, the Lyapunov functional associated with the 
generalized Smoluchowski equation (59) will coincide with the simplified form of J obtained 
by using the Local Thermodynamical Equilibrium condition (53), with u = 0, to express 
J[f] as a functional of p. Using Eq. (55), the energy (3) can be written 

On the other hand, it is shown in [9] that when / is given by Eq. (53), the entropy (4) can 
be rewritten as 

S=^pJ pdh + (3 J \pd 3 r. (66) 
Therefore, the free energy reads 

J = 13 J p d 3 r + 13 J Xpd 3 r P$d 3 r. (67) 

It is shown furthermore in [9] that Eq. (67) can be put in the equivalent form 

J =-/3 [ p f P ^-dp'dh - \(3 [ p<$>d 3 v + C. (68) 
J Jo p 2 J 



This is the Lyapunov functional of the generalized Smoluchowski-Poisson system. Indeed, a 
straightforward calculation yields 

J = I3 (— (Vp + pV$) 2 rf 3 r > 0. (69) 
J KP 

Let us consider particular cases. For the Boltzmann entropy (9), p — pT, and 

J = -Jplnpdh- ^(3jp<5>d 3 r. (70) 

Strictly speaking, J P p(p')/ p' 2 dp' diverges logarithmically as p' — > 0. This means that the 
general formula (68) is only marginally correct for an isothermal equation of state. However, 
assuming a Boltzmann distribution (10) since the begining (see [9]), we can check that 
Eq. (70) is indeed the right formula. Finally, for the Tsallis entropy (13), p = Kp 1 with 
7 = 1 + 1/n, so that 

J=-f3nJpd 3 r -^(3jp$d 3 r + C n , (71) 

where the constant can depend on n. Assuming a Tsallis distribution (14) since the begining 
(see [9]), we find that the free energy reads 
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J = 



7 



-J(PKp^-p)dh- l -pJp§dh. 



This expression is consistent with Eq. (71) and it reduces to Eq. (70) in the limit n 
However, it is more convenient to write the free energy in the form 



(72) 



+oo. 



(73) 



which is also consistent with Eqs. (71) and (70). Under this form, we can say that the high 
friction limit of a g-Tsallis free energy in phase space is a 7-Tsallis free energy in configuration 
space with 7 = (5g — 3)/(3q — 1). We note also that —J//3 can be written "£ — KS" which 
is similar to a free energy with K playing the role of a "polytropic temperature" (see [8]). 
It is not clear whether this formal result bears more physical content than is apparent at 
first sights. 

The generalized Smoluchowski equation (59) can also be written in the form 



dt 



Mp'ipjVp + pV®) 



If we introduce a convex function C(p) through the relation 

pC"(p)=Pp'(p), 

we have equivalently 

where we have defined D = Noting that Eq. (75) is equivalent to 

rp p(p') 



P 



-dp' + dp + C 2 , 



c( P ) = pp r 

Jo 

the Lyapunov functional (68) can be rewritten 

J = -J C{p)d\ - ±0 J p<S>d 3 r, 



(74) 



(75) 



(76) 



(77) 



(78) 



where the term proportional to the mass has been droped. The stationary solution of the 
generalized Smoluchowski equation (76) is given by 



C\p) = - a. 



(79) 



It can be obtained by maximizing the free energy (78) at fixed mass and temperature. 
Similarly, the generalized Smoluchowski equation (76) can be obtained by maximizing the 
rate of free energy production J at fixed mass and temperature (see Sec. IV C). 

The case of Boltzmann and Tsallis entropies in configuration space has been discussed 
previously. The Fermi-Dirac entropy in configuration space S[p] = — J{p In p+(po—p) ln(po — 
p)}d 3 r leads to an equation of state p(p) = — Tm(l — p/po)- For T — > +00, p = pT and for 
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T — > 0, p is a step function with p ~ p in the core and p = outside. The corresponding 
generalized Smoluchowski equation can be written 



dp 
dt 



L>'(Vp + /3p(p -p)V$) 



(80) 



Coming back to the general case and expliciting the relation between <3> and p, we can 
rewrite Eq. (76) in the form 



dt 



(81) 



'd{ P C"{p)V P + Pp J p (r',0^(r-r')dV 
The stationary solution of this equation is determined by the integro-differential equation 

C\p) = -pf p(r>(r - r')d 3 r' - a. (82) 



Equation (81) is one of the most important result of this paper. The general study of this 
equation appears to be of considerable interest in view of its different potential applications. 
This project has been initiated by [28-30] in particular cases. 



C. Generalized stability analysis 

The equivalence between dynamical and thermodynamical stability for the generalized 
Smoluchowski equation (59) is proved in Appendix A. In fact, we can go further and reduce 
the stability problem to an eigenvalue equation as was done in the special case of isothermal 
and polytropic distributions [28,30]. Let p and <3> refer to a stationary solution of Eq. (59) 
and consider a small perturbation Sp that conserves mass. We restrict ourselves to spherically 
symmetric perturbations (non spherically symmetric perturbations do not induce instability 
for non- rotating bodies). Writing 5p ~ e xt and expanding Eq. (59) to first order, we find 
that 



XA 1 d 

r z dr 



r 1 fdSp rf$ d5§\ 



£ V dr ~^ ^ dr ~^ ^ dr ) 
It is convenient to introduce the notation 



(83) 



47rr / dr 

Physically, q represents the mass perturbation q[r) = 5M(r) = J r A^r' 2 8p{r')dr' within 
the sphere of radius r. It satisfies therefore the boundary conditions q(0) = q{R) = 0. 
Substituting Eq. (84) in Eq. (83) and integrating, we obtain 



„2 



d fp'{p)dq\ | 1 dqd® | i7rp d$® ^ 
dr\ r 2 dry r 2 dr dr dr ' 



where we have used q(0) = to eliminate the constant of integration. Using the condition 
of hydrostatic equilibrium Tdp/dr + pd^/dr = and the Gauss theorem in perturbed form 
d5&/dr = Gq/r 2 , we can rewrite Eq. (85) as 
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1 d 



Airpr 2 



( P'{p) dq 



1 1 dq dp Gq 



or, alternatively, 



Anp dr \ r 2 dr J Anp 2 r 2 dr dr r 2 



d ( p'(p) dq^ | Gq = A{ 



(86) 



dr \Airpr 2 dr 



Airpr 2 ' 



(87) 



with q(0) = q(R) = 0. 

Let us now determine the second order variations of the generalized free energy J. From 
Eq. (68), we find that 



S 2 J = -(3 J ^-^pfdh -\(3j 5 P 5*d\ 



(88) 



which must be negative (with respect to mass preserving perturbations) for generalized 
thermodynamical stability in the canonical ensemble. Adapting a procedure similar to that 
followed in [31,29,9], we rewrite Eq. (88) in the form 



(89) 



S 2 J = ( R "dr - -(3 [ R %$dr 

Jo 8irpr 2 \drj 2 Jo dr 



8-Kpr 

Integrating by parts and using the boundary conditions on q, we get 



S 2 J = p q— 



R 



d / p'(p) dq 



lo ^ dr \8ii pr 2 dr 
Using the Gauss theorem, we find 



« d5§ , 
q— — dr. 

dr 



(90) 



S 2j_ p f R q d ( p '( p ) dq )dr | 1 (3 [ R ° q2 dr 
Jo dr\8irpr 2 dr) 2 Jo r 2 



(91) 



or, equivalently, 



5 2 J = -(3 / drq 
2 Jo 



G d / p'(p) d 
r 2 drKAnpr 2 dr 



q. 



(92) 



The second order variations of free energy will be positive (implying instability) if the dif- 
ferential operator which occurs in the integral has positive eigenvalues. We need therefore 
to consider the eigenvalue problem 



d ( p\ P ) d^ | G 
dr \Anpr 2 dr 



q x (r) = \q x (r), 



(93) 



with q\(0) = q\(R) = 0. If all the eigenvalues A are negative, then the critical point is 
a maximum of free energy. If at least one eigenvalue is positive, the critical point is an 
unstable saddle point. The point of marginal stability in the series of equilibria is determined 
by the condition that the largest eigenvalue is equal to zero A = 0. We thus have to solve 
the differential equation 
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±(rt±w\ + GF(ri 

dr \4irpr 2 dr ) r 2 

with F(0) = F(R) = 0. 

We note that the eigenvalue problems (87) and (93) are similar and that they coincide 
for marginal stability in continuity with our previous studies [28-30]. We have also found 
a similar eigenvalue equation by analyzing the stability of barotropic stars with respect to 
the Euler- Jeans equations [31,8,9]. These eigenvalue equations have been solved analytically 
(or by using simple graphical constructions) for an isothermal and a polytropic equation of 
state in [31,8,9,29,30]. It is found that the case of marginal stability (A = 0) coincides with 
the point of minimum generalized temperature 1//3 as predicted by classical turning point 
arguments in the canonical ensemble [18]. The structure of the perturbation profile that 
triggers the instability (in particular the number of nodes) has also been determined in our 
previous papers. The present analysis shows that the structure of the mathematical problem 
remains the same for a general equation of state p = p(p) even if the solutions cannot be 
obtained analytically. 



IV. APPLICATION TO STELLAR DYNAMICS AND 2D TURBULENCE 

A. Violent relaxation and metaequilibrium state 

We consider a two-dimensional incompressible and inviscid flow evolving in a plane per- 
pendicular to the direction z. Let u = — z x Vip denote the velocity field satisfying the 
incompressibility condition V ■ u = 0. The streamfunction ip is related to the vorticity 
wz = V x u by the Poisson equation = — uj. More generally, we can consider a relation 
of the form ip(r) = / g(r — r')uo(r')d 2 r' like, e.g., in the quasi-geostrophic model [32]. We 
assume that the dynamics is governed by the 2D Euler equation 

^ + uVw = 0. (95) 
at 

This equation describes the inviscid evolution of continuous vorticity flows. It also describes 
the mean-field evolution of N 3> 1 point vortices before discrete correlations have developed 
(Vlasov limit) [33,2]. The 2D Euler equation conserves the circulation 

r = J udh, (96) 

the energy 

E=\J u>W\ (97) 

and the vorticity moments (equivalent to the Casimirs) T n = j u n d 2 r. It also conserves the 
angular momentum L — J uor 2 d 2 r in a circular domain and the impulse P — J uyd 2 r in a 
channel (or in an infinite domain). 

When coupled to the Poisson equation, the 2D Euler equation develops very complex 
filaments as a result of a mixing process. In this sense, the fine-grained vorticity u(r, i) never 
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converges towards a stationary solution. However, if we introduce a coarse-graining proce- 
dure, the coarse-grained vorticity u;(r, t) rapidly relaxes towards a metaequilibrium state. 
This is called violent relaxation, chaotic mixing or inviscid relaxation. This "collisionless 
relaxation" is common to all Hamiltonian systems with long range interactions (stellar sys- 
tems, point vortices, non- neutral plasmas, HMF model,...). It physically differs from the 
"collisional relaxation" (or viscous decay) which takes place for much longer times (often 
irrelevant). This implies that the order of the limits t — > +00 and iV — > +00 (or — > 0) is 
not interchangeable [2]. 

There has been some attempts to describe the metaequilibrium state in terms of statis- 
tical mechanics [13,34-36]. In the context of 2D hydrodynamics, the statistical metaequi- 
librium state maximizes the mixing entropy 

S[p] = - J plnp dhda, (98) 

while conserving circulation, energy and all the Casimirs. The mixing entropy is the Boltz- 
mann entropy for p(r, a), the density probability of finding the value u = a in r. The most 
probable distribution is the Gibbs state 

p{*,°) = ^x{°)e-^ +ol) \ (99) 

where x{ a ) = exp(— ^ n >i a n^ n ) accounts for the conservation of the fragile moments r n>1 = 
/ pa n dad 2 r and a,/3 are the usual Lagrange multipliers for T and E (robust integrals) [9]. 
The "partition function" Z — J x(o')e~^ +a ^ da is determined by the local normalization 
condition / pda = 1. The equilibrium coarse-grained vorticity u; = / pa da can be expressed 
as 

57 = = *w + «) = m- (100) 

Since VuJ = f'{ip)Vip and u = — z x Vip, this is a stationary solution of the 2D Euler 
equation. Taking the derivative of Eq. (100), it is easy to show that [37]: 

tj'(^) = -Qui, u 2 = J p{a - ufda > 0, (101) 

where 0U2 is the centered local enstrophy. We note that the relation uJ = f{ip) is always 
monotonic, increasing at negative temperatures and decreasing at positive temperatures. 
Therefore, the coarse-grained vorticity uJ extremizes a if-function [10]: 

S = - J C{uj)d 2 r, (102) 

at fixed circulation and energy, where C{uj) is a convex function, i.e. C"{uj) > 0. Indeed, 
introducing appropriate Lagrange multipliers and writing the variational principle in the 
form 

5S - (35E - a5T = 0, (103) 
we find that the critical points of entropy at fixed circulation and energy are given by 
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C'(u) = - a. 



(104) 



The conservation of angular momentum and impulse can be easily included in the variational 
principle (103) by introducing appropriate Lagrange multipliers Q and U . Equation (104) 
remains valid provided that ip is replaced by the relative streamfunction ip' — ip + ^r 2 — Uy. 
Equation (104) can be written equivalently as 

u = F(Pip + a), (105) 

where F(x) = (C')~ 1 (—x). From the identity 

= -p/C"(u) (106) 

resulting from Eq. (104), uo(ip) is monotonically decreasing if /3 > and monotonically in- 
creasing if (3 < 0. Therefore, for any Gibbs state of the form (100), there exists a H-function 
of the form (102) that the coarse-grained vorticity uJ maximizes (at fixed T, E). It can be 
shown furthermore that ZJ maximizes this functional. We note that C{uS) is a non-universal 
function which depends on the initial conditions. In general, S[uj] is not the ordinary Boltz- 
mann entropy Sb [w] = — J w In ud 2 r due to fine-grained constraints (Casimirs) that modify 
the form of entropy that we would naively expect. In the two-levels approximation {o"o, 0} 
of the statistical theory, S\uJ] is the Fermi-Dirac entropy 

S F . D . [57] = - /(- In ^ + ( 1 - ^ )\n (l - }d\ (107) 
leading to the distribution 



57=- (108) 

In the dilute limit uJ <^ ctq, one recovers the Boltzmann entropy 

S B p} = - J uj\nujd 2 r, (109) 

and the isothermal vortex 

ZJ=Ae~W. (110) 

Other forms of H-functions compatible with the statistical prediction (100) are presented in 
Sec. IV F (from now on, we drop the bar on uj except in case of ambiguity). 



B. Dynamical stability and thermodynamical analogy 

Unfortunately, the statistical theory of violent relaxation is not very predictive because 
the initial conditions are not known in practice and the Casimirs cannot be determined from 
the coarse-grained field once the vorticity has mixed (since a/™ ^ of 1 ) [38]. In addition, the 
high order moments of uj are altered by non-ideal effects (viscosity, forcing,...) so that their 
strict conservation is abusive [39,12]. Finally, the relaxation is in general incomplete so that 
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the ergodic hypothesis which sustains the statistical theory is not fulfilled everywhere [40,7]. 
One aspect of incomplete violent relaxation is that the metaequilibrium structures (vortices, 
galaxies,...) are more confined than predicted by statistical mechanics. 

The only thing that we know for sure is that the metaequilibrium state reached by the 
system is a dynamically stable stationary solution of the 2D Euler-Poisson systems (on a 
coarse-grained scale). Ellis and collaborators [12] have shown that a strong (nonlinear) 
condition of dynamical stability is that u maximizes a H-function at fixed circulation and 
energy. This condition of dynamical stability can be written 

S 2 J=-J C"(u)^^d 2 v - l -f3j 5u5^d 2 r < 0, 

V 5lu I 5E = 5T = 0. (Ill) 

This is similar to a condition of microcanonical stability in thermodynamics. It is therefore 
relevant to develop a thermodynamical analogy and use the same vocabulary as in thermo- 
dynamics to analyze the dynamical stability of 2D flows. In this analogy, S[uj] can be called 
a generalized entropy, (3 a generalized inverse temperature,... In addition, we can introduce 
a Legendre transform J = S — (3E which is similar to a free energy in thermodynamics. The 
condition that u is a maximum of J at fixed temperature and circulation can be written 

C"(cu)^^d 2 r - -13 [ 5cu5^d 2 r < 0, 

V 5u | 5T = 0. (112) 

This is similar to a condition of canonical stability in thermodynamics. The criteria (111) 
and (112) are not equivalent if the "caloric curve" /3(E) presents turning points or bifurca- 
tions [41,12]. This corresponds to a situation of ensemble inequivalence in thermodynamics. 
Since canonical stability implies microcanonical stability (but not the converse), the stability 
criterion (111) is stronger than (112). This has important implications in geophysical and 
jovian fluid dynamics [12,42]. 

The H-function S[u>] maximized by the system at metaequilibrium is non-universal. It 
depends on the initial conditions and on the strength of mixing. The Tsallis entropy 

S q [uj] = --^ [ J(^-uo)d 2 v ) (113) 

is a particular H-function that sometimes occurs in 2D hydrodynamics [4]. However, this 
is essentially fortuitous, as discussed in [7]. In addition, in this context, Tsallis functional 
is not a true entropy. Its maximization at fixed circulation and energy is a condition of 
(nonlinear) dynamical stability not a condition of thermodynamical stability. It leads to a 
particular class of stable stationary solutions of the 2D Euler equation characterized by 

u = A{X-^j) n , (114) 

with A—[(q — l)/3/q]~ and A = [1 — (q — 1)«]/ (q — We shall call this class of vortices 
polytropic vortices by analogy with stellar polytropes [26]. The index n of the polytrope is 
related to the parameter q by the relation n — l/(q — 1). Isothermal vortices are recovered 
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in the limit q — > 1 (i.e. n — > +00). For q = 2, i.e. n = 1, the relationship between u and -0 
is linear [4]. 

To conclude this section, we note that dynamical stability results similar to (111) and 
(112) have been obtained for the Vlasov-Poisson system [11,9]. Therefore, the generalized 
thermodynamical formalism developed in Sec. II can be used to study the (nonlinear) 
dynamical stability of collisionless stellar systems in astrophysics and other systems with 
long-range interactions described by the Vlasov equation. 



C. Relaxation equations for 2D flows 

According to the previous discussion, an important problem in 2D hydrodynamics is 
to construct particular stationary solutions of the 2D Euler equation with strong stability 
properties. We shall consider solutions that maximize a H- function S[u] at fixed circulation 
T and energy E [12]. The explicit construction of such solutions is non-trivial. Expoiting 
the thermodynamical analogy discussed previously and using the MEPP, we can obtain 
relaxation equations (similar to Fokker-Planck equations) that can be used as numerical 
algorithms to construct arbitrary stable stationary solutions of the 2D Euler equation. To 
apply the MEPP, we write the relaxation equation in the form 

duo 

— + u- Vlu = -V- J„, (115) 

where the diffusion current J w has to be determined. The form of Eq. (115) ensures the 
conservation of circulation provided that J w ■ n = on the domain boundary (with normal 
vector n). From Eqs. (97), (102) and (115), it is easy to put the time variations of energy 
and entropy in the form 

E = J J„- V^ 2 r, (116) 

S = - f C"(w)J u • Vcurf 2 r, (117) 

where we have used straightforward integrations by parts. Following the MEPP, we now 
determine the optimal current J w which maximizes the rate of entropy production (117) 
while satisfying the conservation of energy E — and the constraint 

Jl<C(v,t). (118) 

This maximization problem leads to the optimal current 



JL = —D 



ujC"(uu)Vuj + (3(t)ojVi/j 



(119) 



The time evolution of the Lagrange multiplier /3(t) is determined by introducing Eq. (119) 
in the energy constraint (116), using E — 0. This yields 

_ JM"HV,.V^r 
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Introducing the optimal current (119) in Eq. (115), we obtain a relaxation equation of the 
form 



9F + U 



Vuj = VlD 



uC"{uo)Vuj + (3{t)uVi/j 



(121) 



The first term is a generalized diffusion and the second term is a drift. The function (3(t) can 
be considered as a time dependant inverse temperature (possibly negative). It evolves with 
time so as to conserve the total energy E (microcanonical description). The drift coefficient 
£ = Dj3 is a generalized Einstein relation. We shall use Eq. (121) when oo > 0. When oo can 
be positive and negative, we shall prefer the alternative form 



+ u ■ Vcj = V< D' 



(122) 



Equation (122) is obtained from Eq. (121) by setting D' = DujC"{uj). It is straightforward 
to check that Eq. (121) with the constraint (120) satisfies a if-theorem for the generalized 
entropy (102). Indeed, Eq. (117), (116) and (119) lead to 

a=/£ft, (123) 

where we have used E — 0. If uo > 0, then Eq. (123) is positive provided that D > 0. If we 
use the alternative equation (122), we have Duj = D'/C"(iv) so that Eq. (123) is positive 
whatever the sign of ui provided that D' > 0. At equilibrium 5 = 0, hence = 0, which is 
equivalent to 

VC"(w) + (3Vi/j = 0. (124) 

Integrating, we get 

C'(u) = -(3ip - a. (125) 

which returns Eq. (104) with (3 = lim t ^ +00 /3(t). Therefore, as expected, a stationary 
solution of Eq. (121) extremizes the generalized entropy (102) at fixed energy and circulation. 
In addition, it is shown in Appendix B that the relation (34) remains valid so that only 
maxima of S at fixed E and Y are selected by the relaxation equation. The relaxation 
equation appropriate to the canonical situation is obtained by maximizing J = S — j3E with 
the constraint (118). This leads again to an optimal current of the form (119) but with 
constant (3. In addition, Eqs. (123) and (34) remain valid with J in place of S. Hence, the 
free energy J increases monotonically until a state of maximum free energy is reached. 

To conclude this section, we shall discuss particular cases. For the Boltzmann entropy, 
C"(u) = 1/u and Eq. (121) reduces to 

du 



dt 



u ■ Vuo = V 



D 



For the Fermi-Dirac entropy, C"(lo) = l/u(a 
duo 
~dt 



+ u • Vuj = V< D' 



- uo) and Eq. (122 
Voj + (3uj{(Tq — oj)Vip 



yields 



For the Tsallis entropy, C"{oj) = qu q 2 and Eq. (121) becomes 



duj 



D 



(126) 



(127) 



(128) 
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D. A physical numerical algorithm 



Due to the thermodynamical analogy, the relaxation equations proposed in Sec. IV C can 
provide a powerful numerical algorithm to compute arbitrary dynamically stable stationary 
solutions of the 2D Euler equation. Since we are only interested by the stationary solution, 
we can forget the advective term in Eq. (122) and fix D to an arbitrary positive constant. 
We propose therefore the physical numerical algorithm 



(129) 



These equations satisfy the conservation of circulation and energy (robust integrals) and 
increase the generalized entropy (102) until the system has reached a maximum of S at 
fixed T and E. We have seen indeed that a minimum or a saddle point of S[u] are linearly 
unstable via Eqs. (129)-(130). Note that specifying C(uo) does not directly determine the 
equilibrium state because many bifurcations can occur in parameter space (r, E). There can 
also exist local entropy maxima (similar to metastable states in thermodynamics) leading to 
a complicated notion of basin of attraction (see [28] in a related context). These equations 
could be used to study dynamical stability problems on a new angle [12]. Indeed, if a 
stationary solution of Eq. (129) is stable at fixed inverse temperature (3 or at fixed energy 
E, then it is nonlinearly stable with respect to the 2D Euler equation. On the other hand, 
if a stationary solution of Eq. (129) is stable at fixed energy E (variable (3{t)) but not at 
fixed inverse temperature f3, then it violates Arnold's sufficient conditions of stability (see 
Appendix C) but is, however, dynamically stable with respect to the 2D Euler equation. 
The relaxation equations proposed in Sec. II could be used similarly as numerical algorithms 
to construct nonlinearly dynamically stable stationary solutions of the Vlasov equation. 



E. A simplified parametrization of 2D turbulence 

The generalized Fokker-Planck equations derived in Sec. IV C can also provide a simpli- 
fied parametrization of 2D turbulence. The thermodynamical parametrization proposed by 
Robert & Sommeria [43] can be written 



Vp + /3(t)p(a-p)W 



(131) 



where p(r, a, t) denotes the density probability of finding the vorticity level a in r at time t. 
Equation (131) incorporates a turbulent viscosity D and an additional term interpreted as a 
systematic drift [44]. The drift is due to the inhomogeneity of the medium and is supported 
by arguments of kinetic theory in simplified models (point vortices, quasilinear approxima- 
tion,...) [44,33,45]. Usual parametrizations including a single turbulent viscosity correspond 
to the infinite temperature limit {(3 = 0) of the thermodynamical parametrization. Equation 
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(131) increases the mixing entropy S[p] while conserving the energy and all the Casimirs. 
For t — > +00, the solution converges to the Gibbs state (99). 

The equations of Robert & Sommeria [43] are complicated because they take into account 
the conservation of all the Casimirs. This clearly leads to practical difficulties. This also 
leads to physical difficulties because the strict conservation of all the Casimirs is abusive as 
discussed in IV C. We could try to simplify the problem by writing a hierarchy of equations 
for the moments of p. The first equation of this hierarchy is 



9uJ „ „ , T 
— + u • VcJ = V<! -D 

at 



Vu + (3(t)u 2 Vil) 



(132) 



where u 2 is the local centered enstrophy defined in Eq. (101). However, we are now led 
to a difficult closure problem. Kazantsev et al. [46] have proposed to close the hierarchy of 
equations by a Gaussian approximation. This leads to an equilibrium state corresponding to 
a minimum enstrophy state. Although this state may be relevant in some particular oceanic 
problems (Fofonoff flows), this is not expected to be general. 

In this paper, we propose to close the hierarchy of moment equations by a relation of 
the form 



1 



UJ 2 



C"{uj) ' 



(133) 



which can be deduced from Eqs. (101) and (106). This relation is valid at statistical 
equilibrium but we shall use it out of equilibrium as a convenient approximation. In a sense, 
this supposes that the Lagrange multipliers a n associated to the fragile constraints r n> i in 
the Gibbs state (99) can be treated canonically [12]. Substituting Eq. (133) in Eq. (132) 
we obtain the simplified parametrization 



9oj „ „ . r 
— + u • Voj = V\D 

at 



(134) 



which coincides with Eqs. (122) and (120). The function C{uS) is a non-universal function 
which encapsulates the complexity of the fine-grained dynamics and which depends on the 
situation contemplated (ocean dynamics, jovian atmosphere, decaying 2D turbulence...). 
For a given physical situation, we propose to select a form of C{uS) a priori and compute 
the corresponding equilibrium state. Then, we can check a posteriori whether it was a good 
choice by comparing the result with the information that we have on the system. This is 
similar to the notion of prior vorticity distribution proposed in [12]. Of course, this approach 
is essentially phenomenological and explanatory but it allows to deal with complex situations 
which were previously inaccessible. 

To obtain an operational sub-grid scale parametrization of 2D turbulence, it remains for 
one to specify the value of the diffusion coefficient. Heuristic arguments [47,21] or more 
formal kinetic theory [45] suggest that D = Ke 2 uj^ 2 where K is a constant of order unity 
and e is the scale of unresolved fluctuations. This formula is a relatively direct consequence 
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of the general Taylor expression of the turbulent viscosity [48] and it proved to be relevant 
in oceanic modelling [46]. Furthermore, it has been shown in previous works [47,21] that the 
spatial dependance of the diffusion coefficient is important to take into account the problem 
of incomplete relaxation and the formation of self-confined vortices [49]. With the closure 
relation (133), the diffusion coefficient can be expressed in terms of u as 

Ke 2 

D = . (136) 

The same simplifications could be introduced in the parametrization of the gravitational 
Vlasov-Poisson proposed in [21]. 

F. Classification of generalized entropies 

Our parametrization involves a free function C(u). This indetermination is intrinsic to 
the problem of 2D turbulence and not really a flaw of our approach: there is no univer- 
sal entropy S[u]. In order to reduce this indetermination, we shall argue that generalized 
entropies can be regrouped in "classes of equivalence" with the underlying idea that func- 
tional of the same class will produce similar results. Then, for a given physical situation, 
it will be possible to pick a form of C(u) in the corresponding class of equivalence and use 
it in the parametrization (134)-(136). 

We shall now present typical forms of "generalized entropies" S[u] that appeared in the 
literature. Specifically, we shall prescribe analytical forms of the vorticity distribution x(°~) 
in the Gibbs state (99) and compute the corresponding C(u). These analytical expressions 
can be considered as prototypical examples of more realistic distributions. Consider first 
the case where the fine-grained vorticity takes only two values a and G\ > a so that 
x( a ) = Xod(°~ — °~o) + Xid( a ~ 111 that case, uj{jp) is the distribution 

UJ = (To + " ; s , (137) 



which can also be written as a tanh, 



uj = tanh 

2 2 



(138) 



For this distribution, u> — > <7i, <j when — f3ip —> +oo, — oo. Using Eq. (104), we get 

C{u) =plnp+(l-p)ln(l-p), u = pa 1 + (1 - p)a . (139) 

This implies 

C( U ) = ^— ln(^), C"{.) = - 1 -. (140) 

(Ti — <7 \<7i — UJ/ [UJ — <7 X<7i — °°) 

Substituting Eq. (140) in the general parametrization (134), we obtain 



duj 
~dt 



+ u-Vuj = V^D Vuj + P(t)(uj - a )(a 1 - Lu)Vip j, (141) 
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with a diffusion coefficient 



D = Ke 2 ^/(io-a )(a 1 -uj). (142) 

These equations coincide with the two-levels approximation of the thermodynamical 
parametrization (131) . This two- levels approximation has been used in relation with shear 
layer instability [40], vortex merging [43] and in a model of Jupiter's great red spot [50]. 

Turkington [38] has considered the distribution x(<r) = x f° r °"o < & < ci and %(a) = 
otherwise. The choice a n = for n > 1 in Eq. (99) amounts to maximizing the entropy 
(98) at fixed energy and circulation, neglecting the higher order vorticity moments. If we 
assume for simplicity that cr = —A and = +A then u{jp) is the Langevin function 

to = \L[-\(/3i/> + a)], L(x) = coth(x) - -. (143) 

x 

It does not seem possible to write down the function C (u) explicitly. However, the relation- 
ship (143) is qualitatively similar to (138) so we argue that the Langevin-type model falls in 
the same "class of equivalence" as the Fermi-Dirac-type model (i.e., they will produce the 
same type of equilibrium states and bifurcations). 

Consider now a situation in which the fine-grained vorticity can take three values <7i, a 2 
and so that xip) = 3( a ) + Xi^(°" — +X2^(cr — a 2 ). Consider furthermore the dilute limit 
of the statistical theory in which Z ~ 1. In that case, the probability of the non-zero levels 
is given by p t = Ajexp[— (flip + ct)(?i\ and oo = pi<?i + P2&2- This dilute limit corresponds 
to the point vortex model [51]. If we assume furthermore, for simplicity, that a = —A and 
o"! = +A we get the sinh-Poisson relation 

uj = -2A\sinh(A(/ty + a)). (144) 

For this distribution, oj — > ±oo as — > ±oo. The sinh-relationship is observed in the 
late stages of 2D turbulence when the initial condition is a random vorticity field [52] . This 
is because the vortices that form at intermediate times are very intense and isolated as in 
a point vortex gas. Note that the statistical approach based on the conservation of all the 
Casimirs does not work in that case because of viscous effects [39]. Using Eq. (104), we find 
that 

ff M-***-im)- (146) 

Then, 

C(lo) = j sinh- 1 - iv4A 2 A 2 + cu 2 , (146) 

and 

° ^ ~ A V4A 2 A 2 + J' ^ 147 ^ 

Therefore, in the case of 2D decaying turbulence with random initial conditions, the 
parametrization (134) that we propose is 



27 



+ u • Vuo = vlD + p(t)\y/4A 2 \ 2 + uj 2 Vil> j, (148) 
with a diffusion coefficient 

D = Ke 2 ^(AA 2 \ 2 + u 2 ) 1/A . (149) 
If we assumes a Poissonian weight factor x{°~) — exp(— \a\/q) in Eq. (99), we get 

2q 2 (W + a) 

u = ~ i- q ^ + ar (150) 

Pasmanter [53] has noted that this relationship yields results similar to those obtained with 
the sinh-relationship, so they fall in the same "class of equivalence" . 

If we now assume that the local distribution of vorticity is Gaussian so that x{ a ) — e~ a2a 
in Eq. (99), then the cu — ip relationship is linear [35] and can be written 

u) = -n 2 (Pip + a), (151) 

where = w 2 is a constant equal to the centered variance of the vorticity distribution. 
Such a linear relationship sometimes occurs in geophysics [54]. The choice a n = for n > 2 
is equivalent to maximizing the entropy (98) at fixed energy, circulation and enstrophy 
(neglecting the higher order moments). The function C(u) is 

C(u) = ^-u 2 , (152) 

so that the generalized entropy S[u] is proportional to minus the (coarse-grained) enstrophy 
r 2 — J uj 2 d 2 Y. Therefore, Eq. (151) can also be obtained by minimizing the enstrophy at 
fixed energy and circulation [55]. A linear uj — ip relationship is also obtained in the strong 
mixing limit of the statistical theory providing an inviscid justification of the minimum 
enstrophy principle [56]. 

Another choice of probability distribution has been proposed by Ellis et al. [12] in a 
model of jovian atmosphere where the skewness is expected to play a crucial role. They 
relate x{ a ) m t ne general formula (99) to the gamma density z a ^ 1 e^ z (z < 0) and find that 

1 /2 

where Q 2 is equal to the variance of x( a ) an d 2Af2 2 ' is equal to the skewness of x{°~)- 
For this relationship, u — > — oo as —(flip + a) — > I/AQ2 and w ^ -1/A as — j3ip — > +00. 
Therefore, the gamma model (153) is somewhat intermediate between the tanh model (138) 
and the sinh model (144). Using Eq. (104), we find that the corresponding generalized 
entropy is 



CM 



This yields 



uj — — ln(l + \uj) 
A 



(154) 
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UJ 



n 2 (i + Xuy 



C"(u) 



n 2 (i + \uj) 2 ' 



(155) 



For A — * 0, we recover the gaussian model (151) as a special case of the gamma model. The 
parametrization (134) that we propose for this model is 



— + u • Vw = V<! D 

at 



Voo + p(t)n 2 (i + Xufvi) 



with a diffusion coefficient 



D = Ke 2 ttl /2 {1 + Xuj). 



(156) 



(157) 



For A — > (gaussian approximation), we recover the first moment equation of the 
parametrization used by Kazantsev et al. [46] in a barotropic ocean model. 



V. CONCLUSION 

In this paper, we have introduced a new class of relaxation equations associated with a 
generalized thermodynamical framework. The Fokker-Planck and nonlinear Fokker-Planck 
equations corresponding to Boltzmann and Tsallis entropies are recovered as a special case. 
The potential applications of this new class of relaxation equations is considerable and con- 
cerns various domains of physics such as stellar dynamics, 2D turbulence, plasma physics, 
chemotaxis, porous media etc... For example, these equations can serve as numerical algo- 
rithms to compute arbitrary nonlinearly dynamically stable stationary solutions of Vlasov- 
Poisson or 2D Euler-Poisson systems. Physical applications of these equations and numerical 
simulations will be presented in future works. 

We have also clarified the concept of generalized thermodynamics introduced by Tsal- 
lis [3] and applied to self-gravitating systems and two-dimensional vortices by Plastino & 
Plastino [6] and Boghosian [4]. If we consider a collection of stars or point vortices and 
take the long time limit (t — > +oo at fixed N >> 1), the statistical equilibrium state result- 
ing from a "collisional" evolution is correctly described by the ordinary Boltzmann entropy 
Sb [/] = — / / m fd 3 rd 3 v or Sb[oj] — — S oj In ud 2 r but the thermodynamic limit is unusual 
(and does not correspond to N, V — > +oo with N/V fixed) [2]. In the case of stellar systems, 
there is no true equilibrium state but this corresponds to important physical processes (evap- 
oration, gravothermal catastrophe) not to a breakup of thermodynamics [9]. The kinetic 
theory of stars and point vortices is extremely complicated due to the long-range nature of 
the interactions. Kinetic equations can be derived rigorously by using projection operator 
technics but they are non-Markovian and integro-differential [57,33]. The if -theorem for 
the Boltzmann entropy cannot be proved without further approximations. In addition, the 
diffusion coefficient depends on position, velocity and time and this can lead to a confine- 
ment of the structure in physical or phase space. These complicated effects (space and time 
derealizations) can induce, for intermediate collisional times, a deviation with respect to 
an ideal statistical evolution. However, the infinite time limit should be described by the 
Boltzmann distribution. 

On the other hand, for times t <C t Tdax and iV — > +oo, a collection of stars or point 
vortices (or continuous vorticity flows) can achieve a metaequilibrium state as a result of a 
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"collisionless" (or inviscid) violent relaxation. Note that the domain of validity of the colli- 
sionless regime is huge for such systems since the collisional relaxation time t re i ax ~ \£jftD is 
much larger than the dynamical time to [2]. The correct entropy is the Boltzmann entropy 
S[p] = — / p In pdr/d 3 rd 3 v or S[p] = — J p\npdad 2 r for p, the local distribution of phase lev- 
els or vorticity levels. The statistical metaequilibrium state is obtained by maximizing the 
Boltzmann entropy S[p] while conserving energy, circulation and an infinite class of invari- 
ants called the Casimirs. Alternatively, the coarse-grained distribution function or coarse- 
grained vorticity can also be obtained by maximizing a H-function S[f] — — J C(f)d 3 rd 3 v or 
S[uJ] = — J C(p)d 2 r, where C is a convex function, while conserving only mass (circulation) 
and energy (robust constraints). This functional depends on the initial conditions and is 
therefore non-universal. In this context, "generalized entropies" arise due to the presence 
of fine-grained constraints (Casimirs) that modify the macroscopic form of entropy that 
we would naively expect. This is a particularity of continuous Hamiltonian systems and 
this makes the equilibrium state difficult (if not impossible) to predict. A classification of 
generalized entropies in classes of equivalence can however be attempted. 

Furthermore, it can happen that the mixing process during the collisionless relaxation is 
not sufficient to justify the ergodic hypothesis which sustains the statistical theory. In that 
case, the metaequilibrium state results from an incomplete violent relaxation. One possibility 
to take into account incomplete mixing and non-ergodicity is to introduce additional kinetic 
constraints in the statistical approach. This can be done by using relaxation equations 
with a variable diffusivity related to the local fluctuations of the vorticity or distribution 
function [47,21]. Accordingly, the physical picture that emerges is the following: during 
violent relaxation, the system has the tendency to reach the most mixed state described by 
the Gibbs distribution. However, as it approaches equilibrium, the mixing becomes less and 
less efficient and the system settles on a stationary state which is not the most mixed state. 
Then, the evolution is stopped until other effects (collisions, viscosity,...) come into play. 
The state resulting from an incomplete violent relaxation is a (nonlinearly) dynamically 
stable stationary solution of the Vlasov or 2D Euler equation on the coarse-grained scale. A 
strong condition of stability is that it maximizes a if -function at fixed mass (circulation) and 
energy [11,12]. Since this condition of nonlinear dynamical stability is similar to a condition 
of thermodynamical stability, we can use a thermodynamical analogy to study the dynamical 
stability of stellar systems and 2D vortices. In this analogy, the H-function can be regarded 
as a generalized entropy. The H-function maximized by the system at metaequilibrium 
depends on the initial conditions, on the efficiency of mixing and on non-ideal effects. Tsallis 
entropy is just a particular if -function leading to polytropic distributions. It can sometimes 
provide a good fit of the metaequilibrium state in case of incomplete relaxation [4] but this 
is not general [7]. Most galaxies and 2D vortices are not described by Tsallis polytropic 
distribution. A better model is a composite model with an isothermal core where mixing is 
efficient (q = 1) and a polytropic halo (q ^ 1) where relaxation is incomplete [9]. 

In conclusion, generalized entropies arise when hidden constraints are in action: Casimir 
invariants, kinetic constraints preventing mixing, forcing and dissipation, geometrical struc- 
ture of phase space (fractality) etc... We can either work with the Boltzmann entropy and 
try to take into account these additional constraints or keep only the usual constraints (mass 
and energy) and change the form of entropy. The second possibility, while leading to some 
indeterminations, is often more convenient. 
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APPENDIX A: LINEAR STABILITY ANALYSIS OF THE GENERALIZED 

KRAMERS EQUATION 

Let / be a stationary solution of Eq. (28) and 5f a small perturbation around this 
solution. Let us now linearize Eq. (28) around equilibrium and write the time dependance 
of the perturbation in the form 5f ~ e xt . Noting that V6 • U 6 = 0, we get 



A5/ + V 6 (5/U 6 )+V 6 (/5U 6 ) = 



96 J f 
<9v 



(Al) 



where J/ is given by Eq. (25). Multiplying both side of Eq. (Al) by C"(f)8f and integrating 
over phase space, we obtain 



A J C"U)(5ffdhd\ + J C"(f)5fV 6 (5f\J 6 )d 3 r d 3 v 



+ 



C"(f)5fV 6 (f5V 6 )d 3 rd 3 v 

The second term can be rewritten 
1 



C'\f)5f d -^dhd\. 



h = \\ c"(/)v 6 



(8f) 2 U 6 \d 3 vd\ = --J C'"(f)(5f) 2 V 6 f • V 6 d 3 rd\, 
where we have used an integration by parts. From the relation (6), we find 

C"(/)V 6 / = -pu 6± , 



(A2) 



(A3) 



(A4) 



Hence, J 2 = 0. Using <5U 6 = (0, — V5$) and Eq. (A4), the third term in Eq. (A2) can be 
rewritten 



I 3 = -p[ 5fV 6± 5V 6 d 3 rd 3 v = {3j5fv- V5$d 3 rd 3 v. 

After an integration by parts, the fourth term in Eq. (A2) can be written 

d 



-J 



6J r 



dv 



C"(f)5f 



d 3 rd 3 v, 



(A5) 



(A6) 



or, using Eq. (A4), 



-/ 



d 3 rd\. 



(A7) 



Taking the variation of Eq. (25) and using Eq. (A4), one finds that 
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(A8) 



63 f = ~Df 
Hence, 

U = -JS3 r (^ + S^jdPrdPv = - J i^ld 3 rd\. (A9) 



~Df ' J J Df 

where the last equality follows from the conservation of energy (21). Inserting the foregoing 
relations in Eq. (A2), we obtain 

A J C"(f)(5f) 2 d 3 rd 3 v + f3 J Sfv- V5$d 3 rd 3 v = - J ^j^d 3 vd\. (A10) 

We now multiply both sides of Eq. (Al) by <5<£> and integrate over phase space. After 
straightforward integrations by parts, we get 

A J 5f5<5>d 3 rd\ - J Sfv- V5<5>d 3 rd 3 v = 0. (All) 

Combining this relation with Eq. (A10), we find 

A J C'\f){5ffd 3 Yd 3 v + \p J 5f5$d 3 Yd 3 v = - J ^^d 3 rd 3 v. (A12) 

Now, the left hand side is just proportional to the second order variations of the free energy 
5 2 J, see Eq. (16). On the other hand, recalling that J/ = at equilibrium, the second 
variations of the rate of entropy production (31) are given by 



J ^Wv, (A13) 



Df 

and they are clearly positive. Therefore, Eq. (A12) can be rewritten in the simple form 
(34). 



APPENDIX B: LINEAR STABILITY ANALYSIS OF THE GENERALIZED 

SMOLUCHOWSKI EQUATION 

Let w be a stationary solution of Eq. (121) and Sou a small perturbation around this 
solution. Let us now linearize Eq. (121) around equilibrium and write the time dependance 
of the perturbation in the form 5uo ~ e xt . We get 

X5uj + V(5uju) + V(wk) = -V53 w , (Bl) 

where J w is given by Eq. (119). Multiplying both side of Eq. (Bl) by C"{uj)5uj and 
integrating over the whole domain, we obtain 

A J C"{u)){5u) 2 d 2 Y + j C"(Lu)5uX7(5tuu)d 2 r + 

J C"(u)5uV(uj5u)d 2 r = - J C" \u)5ujV3 w d 2 v. (B2) 
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The second term in Eq. (B2) can be rewritten 



1 



c"V)v 



(5cu) 2 u 



d 2 r = ~\\ C"'(uj)(5uj) 2 Vuj ■ ud 2 r, 



(B3) 



where we have used an integration by parts. From the stationary condition (125), we obtain 

C"(lu)Vlu = -(3Vi/j. (B4) 

Since u • V?/> = we conclude that I 2 — 0. Using Eq. (B4), the third term in Eq. (B2) can 
be rewritten 



I 3 = -p[ 5uj5uVipd 2 r = 13 J 5uuV5^d 2 r. 
After an integration by parts, the fourth term can be written 



/ 



5J W - V 



C"(u)8u 



d 2 r : 



or, using Eq. (B4), 



C'"(uj) 



d 2 r. 



Taking the variation of Eq. (119) and using Eq. (B4), we find that 



S3u 



Hence, 



-Du 



C (id) 



(B5) 



(B6) 



(B7) 



(B8) 



h = -JS3 u - (j£ + (3V5^ + 5/3Vv) d2 r = -J + S3 U PV8^ d 2 r, (B9) 

where we have used the conservation of energy (116) to get the last equality. Inserting these 
results in Eq. (B2), we obtain 

A J C'\uj)(5uj) 2 d 2 Y + pj 5ujuV5i)d 2 r = - J + 53^V5i?jd 2 r. (BIO) 

We now multiply both sides of Eq. (Bl) by Sip and integrate over the domain. After 
straightforward integrations by parts, and using Su ■ V5ip = 0, we find 

A J SujSipd 2 v - j Sujvl- V5^jd 2 r = J ' 63 u ■ V5^d 2 

Combining with Eq. (BIO), we get 

A J C"(uj)(5uj) 2 d 2 r + A/3 1 5u5ipd 2 r = -J 

Now, the left hand side is just proportional to the second order variations of the free energy 
5 2 J, see Eq. (111). On the other hand, recalling that 3 W = at equilibrium, the second 
order variations of the rate of entropy production (123) are given by 



-d 2 r. 



(Bll) 



(B12) 



5 2 S 



(SJ U 



Duo 



-d 2 v, 



(B13) 



and they are clearly positive (see discussion after Eq. (123)). Therefore, Eq. (B12) can be 
rewritten in the simple form (34). 
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APPENDIX C: SUFFICIENT CONDITIONS OF STABILITY AND ARNOLD 

THEOREMS 



According to Eqs. (106) and (111), 5 2 J can be rewritten 



The term is bracket is called Arnold pseudo-energy or Arnold invariant. Using an integration 
by parts, we have equivalently 



2{J -u'(ip) 



J (V#) Vrj. (C2) 



First assume that f3 > 0. By Eq. (106), we see that oj\ip) < 0. Therefore, 5 2 J < and the 
system is stable. If (3 > 0, then u'(ip) > and we cannot conclude directly. Let us introduce 
a set of normalized eigenfunctions <f>i(r) such that — A0j = Aj0j with 0j = on the boundary 
of the domain and / <pi<f)jd 2 r = 5ij. Noting that — / 0jA0j<i 2 r = A, and integrating by parts, 
we find that A, = /(V0i) 2 rf 2 r > 0. Then, we decompose 5u> and 5<p on these eignefunctions. 
Using the Poisson equation uj = —Aip, we have 8ui = J2i 8uJi4>i and 5ip = J2i ^y<t>i- Therefore, 

sus^r = J2 (C3) 

We now label the eigenvalues such that < Ai < A 2 < .... Then A« > Ai and we get 

/ 5cu5i>d 2 r < (Sui) 2 = ^ / (5u) 2 d 2 r. (C4) 

J Ai • Ai J 



Therefore, 



5 2 J 



<-\l 



Ai u'(4>) 



{5ujyd 2 v. (C5) 



As a result, if < uj'(ip) < \±, then <5 2 J < and the system is stable. Using Eq. (106), 
this condition can be rewritten (5 > — Aimin C"(u). These sufficient conditions of stabil- 
ity are called Arnold theorems. They imply canonical and microcanonical stability in the 
thermodynamical analogy. Our general formalism shows that the Arnold theorems can have 
applications in many other situations than just fluid dynamics (see, e.g., [9]). 
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